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Abstract 

We have investigated electron dynamics in top gated graphene by measuring the gate admittance 



CZ3 

of a diffusive graphene capacitor in a broad frequency range as a function of carrier density. The 



density of states, conductivity and diffusion constant are deduced from the low frequency gate 
capacitance, its charging time and their ratio. The admittance evolves from an RCdike to a skin- 



effect response at GHz frequency with a crossover given by the Thouless energy. The scattering time 
is found to be independent of energy in the — 200 meV investigated range at room temperature. 

This is consistent with a random mass model for Dirac Fermions. 
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I. INTRODUCTION 



IK 



io|. 



Electron scattering in graphene is still a debated matter (see |l|-|3( and references therein). 
The situation is complicated in this 2-dimensional crystal by the strong influences of the 
substrate, the surface contamination, the effects of static distortions and phonons 
By contrast investigation technics are rather limited; the main diagnosis relies on the tem- 
perature and carrier density (n c ) dependencies of the conductivity. The former is weak for 
graphene on substrates and attributed to phonons 1CH12I] . The later comes from the depen- 
dence of the scattering time r(kp) on wave number which can be varied in a broad range 
due to efficient electrostatic gate doping. We have listed in TABLE-fl] the most predicted 
behaviors for the r(kp) and the corresponding a(n c ) laws 13Ml9||. They can be classified in 
two categories corresponding to a nearly linear density dependence of the conductivity (or 
r ~ kp), or a constant conductivity (r ~ kp ). The former is exemplified by the charge 
impurity mechanism and the latter by the local impurity model. Intermediate situations 
can be explained by an admixture of both mechanisms giving rise to a sublinear density 
dependence of conductivity. Alternatively sublinearity could be accounted for by a single 
mechanism, like the Dirac-mass disorder associated with a random lifting of sublattice de- 
generacy 23j. Indeed, according to Ref.[15(|. a random Dirac-mass mechanism should give 
t = Const.. It translates into a a oc y/n~ c dependence of the bulk conductivity, which should 
not be confused with the G oc ^Jn~ c dependence of a two terminal conductance in ballistic 
graphene 24j, [25 ] . 

Experimentally a linear dependence of conductivity is reported at low temperature with 
a tendency to sublinearity at high density. It is well explained by a Boltzmann conductivity 
c -1 (7ic) = (n c ep + Co) -1 + p s J8(, with a mobility p = evpr jfikp, a saturation at pj 1 and a 
conductivity minimum o"o. The weak effect of dielectric environment 20] and complementary 
measurements of the quantum scattering time 22| suggest an interpretation in terms of 
resonant scattering. 

In the present work we consider the case of top gated graphene at room temperature. 
We investigate scattering using the rf gate admittance. From the in-phase and out-of- 
phase responses we obtain the thermal averages of the density of states p(Ep) and the bulk 
conductivity cr(Ep) as function of energy Ep = hkpVF (vp — 10 6 m/s). We rely on the 
Einstein relation a = e 2 pT> to deduce the diffusion constant T>(Ep) and the scattering time, 



mechanisms 


scattering time 


conductivity 


ite/. 


local impurity 


T ~ \jkp 


a ~ Const 


[13] 


local impurity 


r ~ In kp/kp 


a ~ lnn c 


[14] 


random Dirac-mass 


t ~ Const 


cr~ ^/n^ 


[15] 


charged impurity 


t ~ kp 


cr ~ n c 


[16] 


resonnant scattering 


t ~ kp In (fcp) 


o" ~ n c In n c 


[17] 


ripples 


; (2fl-l) 


cr~nf 


[18] 


acoustic phonons 


r~fe| 


3/2 


[19] 



TABLE I: Main scattering mechanisms proposed for graphene classified according to the carrier 
density dependence : Fermi wave vector dependence of the scattering time r(kp) and carrier density 
dependence of conductivity o~(n c ). For the ripples model, H ~ 1 is the exponent of the correlation 
function for the lattice distortions. The expression for acoustic phonons scattering corresponds to 
a low temperature limit. 

t = 2T>/v F . We shall focus on the hole-doped regime, Ep = —(0-200) meV, where we find 
that r is energy independent. Finally we discuss the possible origin of this observation. 

II. EXPERIMENTAL PRINCIPLES 



In an rf transport experiment, p{E) can be obtained from the quantum capacitance 
cq(Ep) = e 2 J p(E)(—df /dE)dE, where Ep stands now (and below) for the chemical po- 
tential, cq is the thermal average of p(E) and corresponds to the electronic compressibility, 
X = dn/dEp = cq/c 2 where n = J p(E)f(E)dE. The finite compressibility is responsi- 
ble for a chemical contribution AEp in addition to the electrostatic one in the electronic 
charging energy. For a capacitor charge q = —eAn and bias U, AEp = An/x = —Gq/cQ 
adds the electrostatic term, —eq/c geo , in the total energy — eU = —eq/c g . Here c geo and 
c g are the geometrical and total gate capacitance (per unit area) so that one finally has 
Cg(Ep)' 1 = c~} + cq{Ep)~ x . The quantum capacitance term is generally negligible in met- 
als (cq —)• oo) or in back gated conductors (c geo — > 0). 

Using the above definition, one can rewrite the Einstein relation as a = c q T> [2j. cq(Ep) 
can be separated from the constant geometrical capacitance in graphene thanks to its energy 



dependence 26l-l3l|. At finite temperature the theoretical expression is 26] 



2e 2 k B T , . . E F , n 

c ^^T^ H2 + 2coM ^f )] - (1) 

Using nanometers thick gate oxides, c geo ~ 10 fF/fim 2 becomes comparable to the quantum 
capacitance minimum, cq(0) = 4e 2 fcsTln[2]/7r(^,fir) 2 ~ 10 fF//im 2 at neutrality and room 
temperature. 



Top gated graphene offers the possibility to manipulate ac gate currents 



or dynamical 



characterization 32] which is also useful for technological applications 33|, |34(. It has how- 
ever the drawback of a non-linear n c (V g ) relation between gate charge and voltage due to 
the quantum capacitance contribution that complicates the standard analysis. At T = in 
pure graphene one has n c = n gate -n Q [{l + 2n gate /n Q f- b - 1] where n gate = c geo V g /e (see e.g. 
Ref.[2j) and uq = \(c geo hvF) 2 je A — 2 x 10 11 cm -2 (c geo ~ 10 fF//zm 2 ) entails strong devia- 
tions from linear a(n gate ) behavior. In experiment we shall express the measured quantities 
directly as function of the chemical potential Ep using the relation 

E F {V g ) = e(V g - V np ) - / c Q /(c geo + c Q )dV, (2) 

where V np is the gate voltage at the neutrality point. Top gating gives in principle the 
possibility to characterize charged impurity scattering by observing a cutt-off in the Thomas- 
Fermi screening at wave vectors qpF {<1tf ^ 1.5fcp for graphene between SiC>2 and AlOx) 
smaller than a q gate = l/2t ox [35]. However, with t ox ~ 8 nm, q gate ~ 6 x 10 5 cm -1 still 
remains at the lower end of our investigated wave number range kp = 0-30 x 10 5 cm -1 , 
where thermal and impurity effects are prominent. 

The capacitor geometry of our experiment is described in FigJTJ We have measured 
several graphene samples, both in the two-terminal capacitor and three-terminal transistor 
geometries, obtained by the exfoliation method on a thermally oxidized silicon substrate 
with a resistivity p > 20 kQ cm. AFM inspection, performed prior to gate deposition, shows 
a significant roughness of the silicon oxide, which contrasts with the relative smoothness 
of the area covered by graphene. Data presented here refer to two representative capacitor 
samples : samples E9-Zc and C7-F are flakes of dimensions L x W ~ 3 x 1 /im 2 (see Fig{T]) 
and L x W = 2 x 0.6 /im 2 . The use of high resistivity silicon is important to minimize the 
spurious conductance through the substrate. This parasitic contribution can be very high 
at microwave frequencies, and effectively shunt the contribution of a small graphene sample. 



We use palladium for the drain electrode to minimize contact resistance and a gate oxide 
obtained by multistep oxidation of thin aluminum, for a nominal oxide thickness of about 
t ox ~ 8 nm, and finally gold-gate deposition. The distance L a between the drain and the 
gate is reduced to ~ 200 nm, i.e. ~ 6% of the gate length, to minimize the contribution of 
the access resistance R a . With a permittivity k ~ 7 for AlOx we estimate c geo ~ 8 fF//im 2 . 



The device is inserted in a coplanar wave guide used for RF characterization [32 ]. 

The samples are measured in an RF probe station at room-temperature. Bias tees are 
used to control the DC gate-drain voltage V g . The RF scattering parameters are measured 



with a network analyzer and used to calculate the gate-drain admittance Y(ui) |36|. The 
background contribution, which corresponds to a parallel gate-drain capacitive coupling, 
Cq ~ 1.8 IF, is measured in an identical but dummy structure and subtracted. Linear 
response conditions are secured by probing the device well below the thermal noise floor 
with an excitation voltage V r f ~ 1 — 10 mV. We have measured admittance spectra Y{uj) 
up to 10 GHz in a gate voltage range V g = 0-1 V corresponding to average electron and 
hole densities n c = 0-2 x 10 12 cm -2 . For a quantitative analysis, we shall concentrate in the 
following on the hole region where the contact resistance due to chemical hole doping of the 



palladium can be neglected |37 



III. EXPERIMENTAL RESULTS AND ANALYSIS 

Figf2] shows typical admittance spectra of sample E9-Zc for three representative gate 
voltages corresponding to charge neutrality (a), intermediate (b) and large (c) hole concen- 
trations. The purely capacitive response (Re(Y) ~ 0, Im{Y) ~ juC g ) is observed at low 
frequency. Strong deviations from this limit are seen above 1 GHz which are due to finite 
charge relaxation resistance. Solid lines are fits to the data using a 1-dimensional distributed- 
rc model (FigfTJ-b), with linear capacitance and resistance c = c g W and r = (aW)^ 1 . The 
model gives 



Y{uj) = jlocL x 



tanh(jgL) 



(3) 



JQL 

where q = (—jrcu) ' 5 (re = c g /a) is the wave number of the probing rf field in the capacitor. 

We stress that the 1-dimensional character of the rf probing field penetrating the capacitor 
should not be confused with the 2-dimensional electronic diffusion probed by the rf field. 
This ID field distribution is a powerful simplification in the analysis of the capacitor response 



as compared to the more complicated 2-dimensional current distribution in a high frequency 
drain-source conductivity measurement. 
The low frequency development, 

Y{u) ~ juC g + R g (uC g ) 2 + o(u 3 ) , (4) 

gives the equivalent R-C circuit description which corresponds to an homogeneous charge 
distribution with a capacitance, C g = c g LW, and a resistance, R g = a~ l L/?>W . Eq.pj will 
be used for the experimental determination of c g and a. 

The high frequency limit of Eq. ([3]) corresponds to an inhomogeneous charge distribution 
along the capacitor length corresponding to a penetration over a depth, 5 = (2a/c g u) ' 5 <C L, 
from the capacitor drain edge. The asymptotic admittance, 

y = iI_p) v /2W^J (5) 

has a constant phase arg(y) = tt/4 and a a/cJ- dependent modulus. It is markedly different 
from that of an RC circuit with a frequency independent charging resistance, such as a 
contact resistance R c which gives $t(Y) —¥ R~ x and ^s(Y) — > 0. The observation of the 
asymptotic skin effect regime given by Eq.fjSJ) will be taken as a proof that the measured 
conductivity is the bulk contribution equivalent to that measured in a four terminal dc 
experiment. 

The crossing of real and imaginary parts at u c ~ ^a jc g L 2 defines the cutoff frequency 
of the capacitor. It can be expressed as the sum us c = ^a /c geo L 2 + ^V/L 2 . The first term 
is the cutoff of a classical capacitor, the second is a mesoscopic correction due to the finite 
density of state. The correction dominates at charge neutrality and low temperature in top 
gated devices (a/c geo V = e 2 p(EF)/c geo — > 0). One can note that in these conditions, u c 
becomes a rather direct measurement of the diffusion constant, which is a further example 
of a mesoscopic effect where a microscopic property (the Thouless energy KD/ L 2 ) shows up 
in a macroscopic measurement (the cutoff frequency of a capacitor) . 

As shown in Figj2j both regimes are observed in experiment with a cutoff that increases 
with carrier density and levels off at 3 GHz at neutrality (see Figj3]). In particular, the 
full spectrum of Eq.([3]) can be verified at low hole doping demonstrating the prominent 
contribution of diffusion in charge relaxation. In the analysis we have also included the 
effect of the access region as a serial resistance R a . We approximate R a ~ a~ l L a /W ~ 



0.15.R 9 . The inclusion of R a results in a small correction in the investigated frequency 
range, but becomes prominent at very high frequency with an asymptotic limit Y — > R^ 1 
and arg(F) — > for 5 <C L a . In order to quantify the effect of R a we have added in FigJ2] (b) 
(dashed line) the spectrum obtained taking R a — 0. Sample C7-F shows similar spectra 
with, however, a higher cutoff due to larger diffusion constant and smaller gate length. 

From the fits of the full set of admittance spectra we obtain the gate capacitance C g and 
the conductance R~ l as function of gate voltage as displayed in FigJU The capacitance and 
the conductance show a broad minimum at V™ p = 0.67 V (sample E9-Zc) and VV" V = —0.07 V 
(sample C7-F) which we have identified as the charge neutrality point, shifted from zero by 
chemical doping. The solid lines are theoretical expectations for C~ l (V g ) with C geo = Const. 
and the finite temperature expression in Eq.([T]). Deviation from theory at low carrier density 
in FigJU is likely due to inhomogeneity in the chemical doping which is not taken into 
account in Eq.([T]). Indeed it is more pronounced in sample E9-Zc which has a larger doping. 
The fit at large density gives an accurate determination of the geometrical capacitance 
Cgeo = 6.4 ± 0.5 iF/im -2 (sample E9-Zc) and c geo = 5.9 ± 0.5 fF/im -2 (sample C7-F) in 
agreement with the rough estimates from geometry. In the following C~} is subtracted from 
Cg 1 to obtain Cq. R g shows an excess at large electron density (V g ~ 1 V) which can 
be estimated from the difference of the measured resistance to the electron-hole symmetric 
expectation (dotted guide line in FigJl]). We assign this difference to the contact resistance 
due to the formation of a p-n junction in the access region with R c /W < 300 Ohms//mi. 



This value is consistent with theoretical expectation and experiment [38]. We also checked 
that the effect of R c is also seen in the corresponding admittance spectra (not shown) with an 
increase in the series resistance at large V g . For this reason the electron regime is thereafter 
disregarded. The gate conductance i?" 1 depicted in FigJU is similar to the drain source 
conductance measured in graphene transistors of similar size. The order of magnitude of 
the mobility (/x ~ 4500 cm 2 V _1 s _1 at n c ~ 10 12 cm -2 ) is also typical of room temperature 
behavior |7j]. 

Relying on the overall good quantitative agreement of the admittance spectra with the ID 
rf field model and that of the gate voltage dependence of the capacitance, we proceed below 
to the quantitative analysis and deduce in Figj5] the diffusion constant T> from the ratio 
u* = L 2 /3V of the measured quantum capacitance Cq to the charging resistance R g . In 
both cases we observe a linear dependence <j(cq) which corresponds to an energy-independent 



diffusion constant T> ~ 180 cm 2 /s (sample E9-Zc) and D ~ 540 cm 2 /s (sample C7-F) and 
scattering times rikp) = 2D/v F . This corresponds to scattering lengths VpT ~ 40 nm and 



~ 100 nm in agreement with standard estimates [22]. 

Using the experimental cq(V 9 ) and Eq.(J2]) we can plot cq(E f ) and a(E F ) in FigJHJ Small 
deviations from the theoretical estimate for cq in Eq.([T]) are observed in both samples 
which are due to experimental uncertainties and disorder contribution. Uncertainties are 
yet too large for a quantitative estimate of the disorder contribution. The precision of 
the high-frequency compressibility measurement can be increased by using larger flakes 
(more accurate deembedding), thinner gate oxide (smaller electrostatic gate impedance) 
and working at low temperatures (larger contrast in cq). The conductivity data merely 
follow the density of states reflecting the fact that T> = Const. 

IV. DISCUSSION 

We have observed a strong sublinear dependence of conductivity in a broad carrier density 
range at room temperature. Direct comparison between conductivity and compressibility 
in the capacitor geometry at room temperature shows that this sublinear behavior extends 
to the low density limit and is well described by an energy-independent scattering time 
or equivalently a oc \fn~ c . This result deviates from standard low temperature behavior 
(r ~ k F ) explained by charged impurity, resonant scattering models, or standard ripples. 
It could be accounted in a ripple scenario by taking a peculiar value for the correlation 
exponent, 2H = 3/2, which deviates from standard expectation H ~ 1 [2, 



18] . Another 
possible mechanism, less discussed in the literature, would be Dirac-mass disorder proposed 



in Ref. 15J and characterized, according to a numerical calculation, by a oc Ep (r = Const.). 
In conclusion, our RF admittance measurement of graphene capacitors provides new in- 
sight into the energy dependence of the transport scattering time. We report on a new 
behavior in top gated graphene at room temperature characterized by a density indepen- 
dent scattering time. Our experiment gives also access to the electronic compressibility via 
the quantum capacitance in fair agreement with theory. It can be implemented at low tem- 
peratures using a cryogenic probe station and used to characterize the different scattering 
mechanisms in more detail. Our study may also prove useful for the design of high frequency 
graphene transistors. 
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FIG. 1: Graphene capacitor made of a monolayer graphene strip coupled to a metallic gate through 
a thin oxide (thickness t ox ~ 8 nm, permittivity n ~ 7). (a) sketch of the layout, (b) SEM picture of 
sample E9-Zc showing the palladium drain, the gold top gate and the outline of the graphene flake. 
(c) 1-dimensional lumped element description of the graphene capacitor, used for data analysis. It 
includes a distributed resistance r = a~ l /W and capacitance c = c g W where c g = c geo c q /{c geo A-c q ) 
is the gate capacitance and c geo (resp. c q ) its geometrical (resp. quantum) contributions. R a is a 
series access resistance. 
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FIG. 2: Admittance spectra of a diffusive graphene capacitor (sample E9-Zc) for three gate voltages 

V g = (0.68,0.48,0.18) V for the (a), (b) and (c) panels. The signature of skin effect (i.e. Re(Y) ~ 

Im(Y) oc y/f) is seen at neutrality (V g = 0.68 V) in panel (a) due to a larger resistivity. At 

all investigated densities, the admittance spectra are accurately fitted using the 1-dimensional 

strip-line model and a small access resistance R a = 0A5R g . We find respectively (C g ,R g ) = 

(12.5, 4.18); (14.0, 2.64); (16.0, 1.28) (£F,kfi). For comparison we have added in panel (b) the fit 

obtained with R a = (dashed lines). 
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FIG. 3: Gate voltage dependence of the cutoff frequency u c /2tt of sample E9-Zc estimated by the 
crossing of the real and imaginary parts of the admittance spectrum Y(co). The hole region shows 
a trend to saturation at neutrality (V np ~ 0.67 V from which we deduce an (under-) estimate of 
the diffusion constant V(V np ) < L 2 uj c /b ~ 400 cm 2 /s. 
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FIG. 4: Gate capacitance C g and charge relaxation conductance RZ 1 of sample E9-Zc deduced 
from fits of the admittance spectra with Eq.([3|). Labels refer to the admittance spectra in Figj2]for 
hole densities n c = (0,0.65,1.7) x 10 12 cm~ 2 . The solid line is a theoretical fit to the data using 
Eq.([T]) with C geo = 19.2 fF. The dashed line (lower curve) is a guide for the eye representing an 
electron-hole symmetric resistance fitted to the hole values. 
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FIG. 5: Inverse gate capacitance as function of gate resistance. The geometrical contribution l/C geo 
is subtracted to isolate the quantum capacitance term. Data correspond to the hole conduction 
regime where contribution from the contact resistance can be neglected. The observed linear 
dependencies indicate an energy independent diffusion constant. The slopes uj* give T> = L 2 /3uj* ~ 
180 cm 2 s _1 (sample E9-Zc) and V ~ 540 cm 2 s _1 (sample C7-F). 
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FIG. 6: Chemical potential dependence of the quantum capacitance (a) and the conductivity (b) 
deduced form data in FigH] using Eq.Q. 
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